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O ABSTRACT 

<N 

We construct spectral decomposition of ID Fokker - Planck differential operator. 
^ This reveal solution of Cauchy problem. We develop fundamental solution of Cauchy 

5 problem and compare it with one obtained by other means in our former work [5]. 



CN 1. Fokker - Planck ID differential operator 

The Fokker - Planck equation in one dimension is 

dn dn d(vn) d 2 n 

o 



where 

n = n(t, x, v) - density; 
t - time variable; 
x - space coordinate; 
v - velocity; 

a - damping coefficient; 
^ k - diffusion coefficient. 

(1) is evolutionary equation with differential operator 



Q 

U 



d<l> d(v0) d 2 </> 
dx dv dv 2 



where = 0(jc, v). 



Differential operator (2) is called Fokker - Planck ID differential operator. This operator is not self- 
adjoint. Its adjoint operator is: 



* dy/ d^ 3V 

L( ¥ ) = V --av- + k — . 



(3) 



We discuss simple cyclic boundary conditions 

<j){x + a, v) = </>{x, v), - oo < v < oo. (4) 

2. Eigenfunctions and eigenvalues of Fokker - Planck ID differential operator 

Eigenfunctions </> of differential operator L must satisfy with differential equation 

dtp d(v0) d 2 
L(fl = ^ = -v^ + a-^+*_!. (5) 
dx dv dv £ 

and boundary conditions (4). 

We search solution of equation (5) in the form 

* = X{x) V(v). (6) 

The separation of variables provides two equations 

X' = cX, (7) 

kV" + a(vVY - cvV = XV, (8) 

where c = const is a separation constant. 
The solution of (7) and (4) is 

X m = e a ■ (9) 

m 

where m is integer and accordingly c = 2ni — . 

a 

We expect that solution V of (8) satisfies with Fourier transform applicability conditions. The trans- 
formation is 

+oo 

V(v)= J F(co)e imt dco, (10) 

— oo 

and inversely is 

+oo 

F((q)=— \ V{v)e~ i0)v dv. (11) 
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Transformation of equation (8) is 

-kco 2 - acoF' - icF' = AF , (12) 



or 



(aco + ic)F' = -{A + kco )F, (13) 



or finally 



F' A + kco 2 



m 



F aco-2n 

a 



Let us denote 



Take exponent of (17) and get 



(14) 



m 

co*=2n — . (15) 
aa 



Then we have 

F' k(co 2 - a)*) + A + kcol k A + kco* 

= = (C0 + 6)*) . (16) 

F a{to-co*) a a(co — co*) 

Integration of (16) gives 

k r, A + kco 2 . 

lnF = (co + co*Y ln(co-m). (17) 

2a a 

Let us denote 

A + kcol 

n = (18) 

a 



F = e-i (m+0 " )2 (co-co,r. (19) 



Form of expression (19) suggests, that n should be nonnegative integer. Another possibilities we dis- 
cuss later (see end of section 3). Eigenvalues of equation (8) are: 



( 27tm \ 




A = -an - k . 


(20) 


\ aa ) 
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In order to perform inverse transform of (19), we use basic properties of Fourier transformation. First 
of all, we use wellknown fact 



F(6)) = e- ao ' -»V( 



v) = — e 4a. 



(21) 



The second fact is shift theorem 



F(co - a) -> e' av V(v). 



(22) 



Thus we get transformation of the first multiplier in (19): 



-u.M-ni f aV — " v 2 f a V 2,t/h/ " , 

F = e v -> y m0 = 2s- - exp(-/©*v) e 2k =\2ji-\ exp( v) e 2k 

k ) { k I aa 



(23) 



These are eigenfunctions for maximum eigenvalue. Another way to get this solution we show in the 
next section. 

The third fact is differentiation theorem 



ia> F — > 



dV 



(24) 



This implies, that 



F = (a - m)"e~2a (m+e " )2 _> y = — - im 

{ av 



i 

a^a 2mm -— v 

2k — exp( v) e 2k 

k I aa 



(25) 



We use usual convention, that operator acts on functions, staying on the right hand side of operator. 

(25) gives eigenfunctions for the rest of eigenvalues. The aim of following calculations is to express 
them through Hermite polynomials. 

For this purpose we firstly write differentiation operator D in (25) in another way: 



d 2nim 

D = ico* = exp( v) 

av aa 



27tim 

, | exp( v) 

dv ) aa 



(26) 



According to before mentioned convention, this expression means, that we firstly multiply function 

2nim 2nim 

by exp( v), then we differentiate and at last we multiply result by exp( v). Two multipliers 

aa aa 
(the right and the left) have inverse values, therefore we can write arbitrary degree of D in very simple 

form: 
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2mm . 

D = exp( v) 

aa 



d \ 2mm 

— exp( v) 

dv J aa 



(27) 



Therefore we write (25) the following 



a \i 2nim ( d Y 
V = H) In - exp( v) — 



Amm -— v 
exp( v) e 2t 



(28) 



Our next task is to transform exponent in (28) 



Art im -—v 2 

exp( v) e 2k = exp 

aa 



a ( r. Artim 2k 



2k I aa a 



■ exp< 



2fe 



v + 



ATtimk \ ( Azmk \ 
t 2 a ) I «r 2 a J 



•• (29) 



This implies 



V = (-0" exp 



a ( 4-Ttmk \ 



2k y a 2 a J 



a \2 2nim ( d Y 

2 *lJ exp( ^T v) UJ exp 



a ( ATtimk \ 
v + 



2k { a 2 a ) 



(30) 



Hermite polynomials are defined as 



(-1)" d" 
p dz n 



(31) 



where p = e z . 



Add scaling factor and obtain from (31) expression for derivative: 



dz n 



e -az = (-iy a n/2 e -az H n (^ Z ) 



(32) 



a 



Substitute a = — and get the following expression 



V = (0"(-)" /2 exp 



a ( Axmk \ 
2k { a 2 a ) 



a \i 2nim 

2n — exp( v) exp 

k I aa 



a f ATtimk \ 
2k { V+ ^x 2 a r ) 



( \~a~( ATtimk Yt 

"•W* v+ ^r (33) 



Combine two exponents, drop constant multipliers and get final expression for eigenfunctions of dif- 
ferential operator (8): 



2nim ( a \ ( \~a~( ATtimk\\ 

: exp(- — v) exp - - - [v + — ]]. (34) 
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Eigenfunctions of differential operator (5) are equal to 



m v ( a t\ ( I a ( Ajtimk\\ 
*- = exp(2^/ - (x - -)) exp [- - v^„^ - + — jj. (35) 



3. Eigenfunctions and eigenvalues of adjoint differential operator 

Eigenfunctions yt of adjoint differential operator L* must satisfy with differential equation (see (3)) 

* dy/ dy/ d 2 y/ 

L (y/) = Ay/ = v av — + k ^— - . (36) 

dx dv dv z 

We perform substitution (7) once again. This time we use definition 

Y m = e- 2 * im -a; (37) 

instead of (9). The difference is in the m sign. 

After variables separation (analogous to (6-7)) we have the following equation (analogous to (8)): 

in 

kW" - avW - 2m — vW = AW. (38) 
a 

This operator is adjoint to (8). Note, that m has another meaning, due to sign changing in (37). This 
fact is crucial for orthogonality analysis (see below). 



Let us perform substitution 



a 



W = exp|^ v 2 ] V(v); (39) 



a 2 i '-- a 



W" = exp| — v 2 V" + 2-vV' + (-v) 2 V+-V 
\2k ) \ k k k 



in (38). Thus V satisfies with (8). Therefore eigenfunctions of adjoint operator (38) are 



( 27tim \ ( I a 
W ffl „=exp - — v \H H Uj- 



v+ ^rJ) (40) 



Substitution (39) does not involve A. Therefore eigenvalues of differential operator (5) and adjoint 
differential operator (37) are the same ( see (20)) in accordance with theory [1,2]. 
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Eigenfunction with maximum eigenvalue for given period a is equal to 



2mm 

W mQ = exp( v) . (41) 

aa 



This solution we obtain in another way (not depending from Fourier transform), equating coefficients 
by degrees of v in (38) to zero: 

aW' mQ + 2m-W m0 = 0; (42) 
a 

kWZ = XW mQ . 

These equations are compatible. The first equation gives solution (41) and the second gives expres- 
sion (20) (with n = 0) for eigenvalue. 

We get corresponding solution of (8) from (42), performing inverse substitution (39). Obviously, it 
coincides with (23). 

Calculate commutator of two operators: operator D from (26) and operator 

d f ( d 



L - k\ — I + a\ — |v - cv - X. We get easily 



L(D(V))-D(L(V)) = -aD(V) (43) 



Therefore L(V) = implies (L + a)(D(V)) = 0. 

This gives the rule for building of L eigenfunction with eigenvalue X - a strarting from known eigen- 
function with eigenvalue X. 

Starting from basic solution (23) we get (28) by induction. This is another way to get expressions for 
eigenfunctions, independent from Fourier transform. 

d Til (X 

Another important differential operator is D_ = 1- 2m 1 — v. This operator inversely builds 

dv aa k 

eigenfunction of L with eigenvalue X + a from known eigenfunction with eigenvalue X. We see, that for 
n = must be D_(V mQ ) = and we can find V m0 from this equation. 

We can represent L as product of D and D_ , analogous to Schroedinger factorization method. All 
theory is similar to the quantum mechanics linear oscillator theory. 



Eigenfunctions of (36) are 



( III V \ ( 




A mink Yj 




y/ mn = exp -2m —(x + —)\HJ-\ 


m 




(44) 






a 1 a J J 





We built all eigenfunctions of ID Fokker - Planck differential operator. We can return to open ques- 
tion - which values can posess n in expression for eigenvalues (20). In order to answer to this question we 
consider equation for eigenfunctions (38) of adjoint differential operator for m = 0. The point v = is regu- 
lar point of this equation, so we can search its solution in the form of Taylor series by degrees of v. We see, 
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that result is 0(e k ) unless X = -na. If X = -na the series breaks and solution is Hermite polynomial. 
4. Orthogonality of Eigenfunctions 



It is wellknown, that eigenfunction of differential operator is orthogonal to eigenfunctions of adjoint 
differential operator, which belong to another eigenvalues. We check below this orthogonality for our spe- 
cial case. 

Orthogonality of V mp and W mq follows from orthogonality of Hermite polynomials. The product 
V W is equal to 



\nim ( a n \ ( \~a~( 



Animk Y, 
+ H, 



a 2 a ) 



v+ — — 11= (45) 



a 2 a )) 





a 


r Anmk ^ 2 




a 


exp 


~2k 


v a 2 a J 


exp 


~2k 










exp 



v + 



Animk \ 

a 2 a J 



( [~a~f Ammk\\ ( [~cT( \nimk\ 

"{ilk v+ ^r Hiikh^a 



a ( \nmk \ 



2k y a 2 a ) 



exp(-£ 2 )//„(£) H q (0 



ra~f \nimk \ 

where c — \l — v H 

\ 2k { a 



We know (see [3], Theorem 54), that 



I pq = J pH p (4)H q (t)dt = S pq 2 p p\^. 



(46) 



We apply parallel translation of integration path on complex plane and get: 



V mp W mq dv = e\p 



( Ajtmk \ 



2k I a 2 a ) 



= exp 



a ( \nmk "f 
2k{ a 2 a J 



[2M 



exp(-Z 2 )H p (t) H q (4)d{ = 



S pq 2 p P \. 



(47) 



(47) proves, that eigenfunctions V mp and W mq , which belong to the same m, are orthogonal. But V pt 
and W qm functions are not orthogonal, when p * q. Multipliers X p and Y ' (see (7) and (37)) are responsible 
for orthogonality in this case. 



X p Y q dx = aS pq . 



(48) 
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This proves orthogonality of </> mp and y/ nq 



\dx\ 



Wnqdv = a5 mn 8 exp 



a ( \nrnk \ 
2k y a 2 a J 



Ink 



a 



2 p p\. 



(49) 



We can write for arbitrary function G(x, v) the following expansion 

m-oo p=oo 

G(x, v) = £ £ A mp <V (50) 

m=-oo /?=0 



where coefficients A mn are equal to 



^Itin — eX P 



a ( \izmk \ 
2k { a 2 a J 



I 1 " °° 

V^^iW G(x ' v) Wmpdv - (51) 



5. Solution of Cauchy problem. Fundamental solution. 

According to previous results, solution of Cauchy problem with initial condition n = G(x, v) is 



m-oo p-oo 

n(t, x,v)= £ Z ex P 

m=-oo p=0 



-(ap + £ 



f 2;rm 
I aa 



)t 



l mprmp-> 



(52) 



where coefficients A mp are given by (51). 



If G(x, v) = S(x - %)8(v - u) then A mp = — exp 



a ( A-Ttmk ^ 
2k y a 2 a ) 



a/^^^,^ and finally 



8{ 



1 r~a~ m ^±°° 

x-^)S(v-u)= -A/— — 2, ex P 
a \ 2nk m= _„ 



or f \nmk \ 



2k y a 2 a ) 



h 2Pp\ 



(53) 



or 



1 / a m=+a = ( {x- \ ( (v - u) \ 

5(x - %)8(v ~u)-—\ Y, ex P 2/n/n exp 2nim x (54) 

a y 2nk y a J y aa J 



p=oo I 

x Y exp 



a ( \izimk\ 
2k y a 2 a J 



( \~a~( Ajtimk\\ ( \~oT( Animk 



In this way we obtain finally expansion of fundamental solution of Fokker - Planck equation: 
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1 I a m=+ °° 
G a (t, x, £,v,u) = — a / — — £ exp 
a y 2nk m= _„ 



( 27tm 
-k 1 1 



y aa 



expl^^Wf^^V (55) 
a J I aa I 



exp(-«/?f) 



or f \izimk \ 
2k { a 2 a J 



H HTk 



a ( Ammk\, 
v+— — II H, 



a ( \nimk Y\ 
u + 



where E, - initial point of space; 
u - initial velocity. 

Index a below points to interval's length. 

We can simplify (55) by explicit summation of the second sum. As an instrument for this task we use 
following identity found in [3] (Theorem 53) with reference to the original work [4]. 



n=oo g 2 

E 



n=0 



2 n n\ 



s n H n {x) H n {y) 



. exp 



2 2/ n2 

x - y (x - ys) 
l-s 2 



(56) 



This identity is true when Lvl < 1. 

We perform substitution in (55) 



exp(-orf) = s; 

I a ( A7timk \ 
i2k{ V+ ^a-) = ^ 

f~a~ ( Ajrimk \ 

i2k{ u+ ^ 



2 a ) 



= v. 



(57) 



After that the second sum in (55) is equal to 



p=oo f P 



E ^—exp[- M 2 ]H p (fi) H p (v) = 
P =o z ' P- 



(58) 



exp 



-fi +v 



Vl-s 2 



exp 



fi 2 -v 2 (fi- vs) 2 
2 l-s 2 



Vl-5 2 



exp 



l-s 2 



or 



1 j a m =+°° 
G a {t, x, |, v, u) = - A/ — — £ exp 
a y 2izk m= _„ 



( 2itm 

-k I / 

^ aa 



exp 2nim 



(*-£) 



-2a 1 



(59) 



. (y-u) 

x exp 2k im lexp 



aa 



(ju - exp(-arf)v) 
1 - exp(-2af) 
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We return from /u, v to v, u variables 



1 l~a~ m ^°° 

G a (t, X, £ V, K) = - A/ — — X ex P 

a y Ink m= _oo 



-k\ f 

y aa J 



exp \2mm 



a j ( (v — uY. 

= exp 2mm x (60) 

^\- e -2at y aa ) 



xexp 



a 



2k(\-e- 2a <) 



(v-e- a, u) + (l-e- c ") 



<\7timk 



Expand square in the last exponent 



1 I a m=+ °° 
G a (t, x, 4, v, u) = - A/ — £ exp 

a i LK K m——oo 



( 2kwi \ 

-k\ \t 

y aa J 



exp Vlnim 



( *-£) 

a ) f_ . iv-u)' 

exp \2nim x (61) 



1 _ e -2«r ^ era 



xexp 



a 

2k(l - e~ 2aV ) 



(v - e~ at uf + 2(v - e-°"u)i\ - e al ) ^= - (1 - £>""') 



A-Ttimk ^ _ atl (\Kmk^\ 



y a 2 a j 



and then arrange arguments of exponents. This gives the final form of fundamental solution of initial value 
problem with simple cyclic boundary conditions (4): 



G a it, x, v, u) = 



x £ ex P 



. 2kwi \ ( 2 1 — e 
-k t 



a [2xk(l - e~ 2at ) 

2nim 



exp 



a (v-e~ a 'u) 2 
2k (l-e- 2 "0 



(62) 



aa { a 1+ e~°" 



exp 



1 2 (v-e- a ' U y 

< X -4) + -(y-u)-- \ J 

a \ a a (1 +e °") 



6. Fundamental solution for the whole space 



Direct interval length a to infinity and get from (62) fundamental solution for the whole space. The 

27tm 

Fourier series sum is transformed to Fourier integral. Formally, we perform substitutions a = and 



dco ~ 



2k 



G co (f, x, v, u) ■ 



1 



2k \2Kk(\ - e- 2a <) 



exp 



a iv-e~ at uf 
2k (l-e-2«0 



(63) 



+oo 

x J" exp 



a { a 1+ e-°" 



exp 



1 2 (v-e~ at u) 

im\ (x - £) + - (v - k) - - 

a a (1 + e at ) 



dm. 



We perform calculations (63) using Gaussian integral 
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1 +<x 

— f 
In J 



e - Dm \ imp da = 



1 -P- 



2AkD 



e 4D. 



(64) 



For our task D and p are equal to 



1 2 (v-e" af w) 

» = (x - c) H (V — It) 



(65) 



* ( 2 l-e-"' 

D= —\t 

a 2 { a 1+ e~ at 



(66) 



Use another variable D*: 



D* 



at {\-e- 2a, )-2(\-e- at f 
2a 4 



(67) 



This variable was introduced in our former work [5] as determinant of quadratic form matrix. We have: 



D = D* 



2ka 
(1 -e- 2at ) ' 



(68) 



Perform transformation and get another form of fundamental solution of Cauchy problem for the 
whole space: 



G co (f , x, |, v, u) ■ 



2nk(\ - e~ 2at ) 



exp 



a (v-e~ at u) 2 
2~k (l-e- 2a <) 



(69) 



1 f(l-e- 2c ") 



2 { 27tkaD* 



1 (l-e- 2c ") 



expi- 



4kD* 2a 



1 2 (y- e -°' M ) 

(x - £) + - (v - u) - - — 

a a (1 + e at ) 



Join two exponents 



G co (f, x,£, v,w) = 



2^(l-e- 2 «0 



(70) 



i f(i- e - 2c ") y 



2 { 2nkaD* 



exp< 



4kD* 



(l-e- 2a <) 
2a 



. 1 , 2 (v f «f (l-e- 2 '")-2(l-e- < ») 2 (v- e ^ M ) 2 

(* - |) + - (V - U) - - + 5 -r— 

a a (l+e at ) J a 5 (1 - e Zal ) 



or 
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G„(/, x, t, v, li) = = x 



(71) 



xexp< 



4kD* 



(l-e- 2at ) 
2a 



(x-Z) + -(v-u) 
a 



- \(x-4)+-(v-u) 

ay a J (1 +e at ) 



+ — (y-e a, u) 



at,,\2 



a L 



In order to perform final step, we introduce another set of variables: 



x = x-(£+-(l-e- a ')); 
a 



(72) 



v = v - ue 



These variables were introduced in [5] as integrals of characteristic equations. Subexpression of (71) 
is equal to 



1 v 

(x - £) H (V - U) — X -\ , 

a a 



(73) 



and whole (71) is now 



G co (f , x, v, m) = == exp • - 



\nk<D* 



4kD* 



(l- e - 2a ') 



2a 



vf 4f v 

x+ — U + -, 

a J a I a )(\+e~ at ) 

v ' y 



(74) 



We see, that exponent's argument is quadratic form of new variables x, v. 



Very simple calculation brings this expression to final form: 



1 1 1 

G co (f , x, |, v, w) = exp ^ - — — 

Ank-4Di 4-kD* 



2a 



(75) 



(l-e- a ')-^(l-e- 2c ") I xv + l -.-^(l-O+^rd-O I v 



f 2 



2«f\ ,~2 



2a 3 



As expected, this result coincide with obtained in [5]. 

We checked our considerations, comparing result with earlier one obtained by other means. 
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DISCUSSION 

The main result of these considerations is the closed form expression for eigenfunctions and eigenvalues of 
ID Fokker - Planck differential operator. The spectral properties of ID Fokker - Planck differential opera- 
tor are very similar to properties of Hermitian quantum mechanics operators, especially - harmonic oscilla- 
tor. All eigenvalues are real. For cyclic boundary conditions the spectrum is discrete, for the whole space 
the spectrum is continuous. Eigenvalues are nonpositive. The maximum (zero) eigenvalue corresponds to 
space uniform Maxwellian velocity distribution with zero mean velocity. Fundamental solution of Cauchy 
problem is obtained and compared with obtained in former work [5]. These results enable solution of more 
complex problems. 
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